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Time-Resolved Photometry of GK Per during the 1996 Outburst 
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Abstract 

GK Per is a unique cataclysmic variable star which has showed a nova explosion as well as dwarf nova- 
type outbursts, and has the intermediate-polar nature. We carried out y-band time-resolved photometry 
and B-band monitoring during the 1996 outburst. This outburst lasted about 60 d and is divided into 
three parts: the slow rise branch for 35 d, the gradual decay branch with a decay rate of 20.0 d mag~^ 
for ~16 d, and the rapid decline branch with a rate of 5.6 d mag~^ for '--^lO d. The B — V color became 
bluest {B — V ^ 0.18) about 10 d before the outburst maximum, which supports an idea that the outburst 
in GK Per is of the inside-out type. The spin pulse, 440-s quasi periodic oscillations (QPOs), and '^5,000-s 
QPOs were detected in our light curve, as previously seen in X-ray and optical observations. In addition, 
we report the discovery of ^^SOO-s periodicity, which is shorter than the spin period. 

Key words: accretion, accretion disks — stars: novae, cataclysmic variables — stars: intermediate 
polar — stars: individual (GK Persei) 



1. Introduction 

Cataclysmic variable stars (CVs) are a group of binary 
stars consisting of a white dwarf (primary star) and a late- 
type secondary star (for a thorough review, Warner 1995). 
Surface gas of the secondary pouring out from its Roche 
lobe is transferred via the inner Lagrange point to the 
primary. While the gas is accreted by the white dwarf 
through an accretion disk in usual CVs, the inner region 
of the disk is truncated or the accretion disk is completely 
obstructed to be formed in CVs containing a white dwarf 
with strong magnetic fields. In this case, the gas is forced 
to move along the magnetic fields and is accreted into 
the magnetic poles on the white dwarf. These magnetic 
CVs are classified into two subgroups of polars having a 
white dwarf with a spin period equal to the orbital period 
and intermediate polars (IPs) containing a white dwarf 
asynchronously rotating (for a basic review, e.g. chapter 
8 and 9 in Hellier 2001). 

CVs show two kinds of large-amplitude sudden bright- 
enings, namely, nova eruptions and dwarf nova- type out- 
bursts. The former are driven by thermonuclear runaway 
reaction of the material accreted on the white dwarf sur- 
face (Starrfield 1999, and references therein). The ampli- 
tude is 8-15 mag. The recurrence cycle is believed to be 
typically 10** yr or more, but this cycle highly depends 
on the primary mass, the mass transfer rate, the chemi- 
cal abundance, and so on. For example, the shortest re- 



currence period of nova eruption in U Sco is 8 yr, with 
recorded outbursts in 1863, 1906, 1936, 1979, 1987 and 
1999. On the contrary, the thermal instability in the ac- 
cretion disk is generally accepted as the driving mecha- 
nism for the dwarf nova-type outburst (for a review of the 
disk instability model, e.g. Osaki 1996). The amplitude 
is 2-8 mag, and the recurrence cycle is several days-a few 
years. 

GK Per is the only star which has the IP nature and 
has shown a nova eruption and dwarf nova-type outbursts. 
On 1996 February 26, a new outburst was caught in this 
enigmatic CV by J. McKenna (vsnet-obs 2211). Following 
his notice, we started V-band time-resolved photometry 
and accompanying B-band monitoring. The results are 
reported in this paper. We briefly summarize the history 
of the GK Per study in the literature in section 2. Section 
3 and 4 are devoted to describe the observation and the 
result. Discussion is made in section 4, and we summarize 
the paper in section 5. 

2. History of the GK Per study 

GK Per is a quite unique star and has cast challenges 
to researchers in various fields. As this star has a long 
research history, we here summarize a few decades. 

GK Per was originally discovered as a nova in 1901 
(see Hale 1901; Pickering 1901). In the decline phase 
of the nova eruption, this nova showed peculiar oscilla- 
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Table 1. Loe of the observations. 
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* Calculated using the ephemeris given by Morales-Rueda et al. (2002): Tq{HJD) = 2450022.3465- 
t The comparison star is HD 21588 {V = 9.069, B - ^=0.210). 
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tions of ~1 mag with a typical cycle of several days. The 
nature of these oscillations is still unknown (Bianchini 
et al. 1992). The distance was derived to be 460 pc 
by study of the gaseous shell ejected by the nova out- 
burst (Warner 1976), although Duerbeck (1981) quoted 
525 pc. Kraft (1964) first found that GK Per is a double- 
line spectroscopic binary with the secondary star of K2 
IVp, showing spectra having a blue continuum with the 
absorption lines of the secondary star and emission lines 
of H, He I, He ii, [O ii], and so on. They deduced the 
orbital eccentricity and the orbital period to be 0.4 and 
1.904 d, respectively. These high eccentricity and the or- 
bital period were supported by Bianchini et al. (1981) and 
Bianchini et al. (1981). Crampton et al. (1986), how- 
ever, showed that a circular orbit fitted the radial veloc- 
ity variation and measured the orbital period of 1.996803 



(±0.000007) d, the primary mass {Mi) of 0.9 Mq, the 
secondary mass (M2) of 0.25 Mq, and the inclination 
< 73°. This long orbital period has been confirmed by 
Reinsch (1994) and Morales-Rueda et al. (2002) who gave 
a new ephemeris used in this paper. The spectral type 
of the secondary has been determined to be KO V-III 
(Garhck et al. 1994), K2V-K3V (Reinsch 1994), or Kl IV 
(Morales-Rueda et al. 2002). Morales-Rueda et al. (2002) 
constrained the lower limits on the mass: Mi > 0.87±0.24 
M0and M2 > 0.48±0.32Mq. Mi was also estimated 
to be O.36-O.86M0, using an X-ray observation by the 
ASCA satellite (Ezuka, Ishida 1999), although Cropper 
et al. (1998) tried to estimate Mi using Ginga data ob- 
tained in 1987, but failed. 

Very Large Array (VLA) observations of GK Per car- 
ried out by Reynolds, Chevalier (1984) revealed non- 
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Fig. 1. (upper panel) The long term light curve of the 1996 
outburst of GK Per. The open circles and the vertical bars 
represent the daily averaged visual magnitudes reported to 
VSNET and the standard deviation. The large values of the 
standard deviation mean large- amplitude modulations of GK 
Per (see text). The small dots are our observations, (lower 
panel) The B — V color. GK Per became bluer during the out- 
burst, and got to the peak about 10 days before the outburst 
maximum. 

thermal radio emission with a partial shell morphology, 
which coincided with the brightest portion of the opti- 
cal shell. Extended emission around GK Per detected in 
the far infrared IRAS imaging was interpreted to be an 
ancient planetary nebula ejected from the central binary 
system (Bode et al. 1987). This was supported by radio, 
far-IR, and optical observations carried out by Seaquist 
et al. (1989), and the 230 GHz observation in the J = 2 ^ 1 
line of ^^CO also supports the suggestion that GK Per is 
surrounded by a /ossiZ planetary nebula (Scott et al. 1994). 
Hessman (1989), however, cast doubt on this interpreta- 
tion by a large-scale survey of the ^^CO dust emission 
around GK Per. The lUE ultraviolet observation made by 
Evans et al. (1992) gave some evidence of abundance gra- 
dients in the nebula. Based on CCD imaging, Anupama, 
Prabhu (1993) suggested that the nova shell is clumpy 
and asymmetric, and the expansion rate is 0.31 (±0.07) 
arcsec yr~^. Dougherty et al. (1996) detected, from IRAS 
60 fim and 100 /xm observations, a far IR emission region 
extended up to 17 arcmin, which is suggestive of a bipolar 
outflow. 



King et al. (1979) identified GK Per with the X-ray 
source A03274-43. This star is a very hard X-ray emitter 
among classical novae (Becker, Marshall 1981). Cordova 
et al. (1981) noted that there was some evidence for vari- 
ability on a 50-s time scale in the Einstein X-ray light 
curve. X-ray variabilities on timescales of 100 s and hours 
were also suggested by Cordova, Mason (1984). Watson 
et al. (1985) finally detected strong, coherent, hard X-ray 
(>2 keV) modulations of a period of 351 s, using EXOSAT 
during an optical outburst in 1983, which established the 
IP classification of GK Per, interpreting the 351-s period 
as the spin period (Pgpin) of the white dwarf. They also 
found that the hard X-ray emission varies aperiodically on 
a timescale of approximately 3000 s. The 351-s pulse was 
observed in quiescence with EXOSAT (Norton et al. 1988) 
and Ginga (Ishida et al. 1992). The pulse profile in quies- 
cence, however, is in a double-peaked shape, different from 
that in outburst. Ultra-high-energy radiation (>10 TeV) 
was not measured from GK Per (Alexandreas et al. 1991). 
Helher, Livio (1994) re-analyzed the EXOSAT data previ- 
ously published by Watson et al. (1985) and found quasi 
periodic oscillations (QPOs) of a timescale of '^5000 s. 
They proposed the disk-overflow accretion model for these 
long-period QPOs. The 351-s periodic modulations and 
aperiodic intensity variations with a typical timescale of 
several thousands of seconds were detected again during 
the 1996 outburst by ASCA (Ishida et al. 1996). 

Regarding optical observations of short term variability, 
Robinson, Nather (1977) reported an initial result that, 
among ten old novae including GK Per, coherent oscilla- 
tions were not detected except in DQ her. QPOs with a 
timescale of ~380 s were subsequently found in quiescence 
(Patterson 1981)^. U band time-resolved photometry by 
Mazeh et al. (1985) gave evidence of '^400-s periodicity in 
two nights during the 1983 outburst, ~360-s one during 
the decline phase near quiescence, and 0.8-h one which 
might be an optical counterpart of the ~3000-s pulses in 
X-ray (Watson et al. 1985). They suggested relationship 
of a beat phenomenon among Pspin, ~400 s, and 0.8 h. 
Patterson (1991) reported the discovery of the spin pe- 
riodicity in U band and QPOs of ^^355 s in quiescence. 
During the 1986 outburst, U-band photometry performed 
by Pezzuto et al. (1996) detected periodicity of 40 s, ^^370 
s, 20-80 min. In the case of spectroscopic observations, 
Reinsch (1994) gave first firm evidence for the variation 
of the V/R ratio of Balmer and He i, ii lines with periods 
of Pspin, 430-500 s, and 2,000-7,000 s, and the equiva- 
lent width (EW) of these lines modulates with a period 
of 2,000-7,000 s (see also Hutchings, Cote 1986). GarHck 
et al. (1994) also found Pspin modulations of the V/R ra- 
tio and EW of Balmer lines. Morales-Rueda et al. (1996) 
drew power maps of frequency versus line velocity of H/? 
and He ii 4686 from spectra acquired during the 1996 out- 
burst, and detected ~5,000-s QPOs. Re-analysis of the 
same data by Moralcs-Rueda et al. (1999) showed that 
Uf] and He ii have a blue-shift bias for the 4,000-6,000-s 

^ Patterson (1981) listed the QPO in a table with reference to 
Patterson, J., 1981, ApJ, in press. This paper, however, have 
not yet been published to the authors best knowledge. 
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Fig. 2. Daily light curves. Various-amplitude modulations, such as QPOs and flickerings, with various timescales were caught. The 
tick marks and the upper arrows point at the calculated maxima of the 0.05612-d and the 0.05947-d oscillations (see text). 
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QPOs in the power spectra of hne flux in discrete veloc- 
ity bins, which supports the disk-curtain beat mechanism 
rather than the disk-overflow accretion model proposed 
referred to above. 

Dwarf nova-type outbursts have been seen in GK Per 
every few years since about 60 years after the nova explo- 
sion. The quiescence magnitude is around my=13.2, but 
has been modulating with a cycle of a couple of decades. 
This may be due to the solar-type activity in the sec- 
ondary star (Ak et al. 2001; Bianchini 1990; Warner 1988). 
These outbursts are characterized by peculiar features, 
a long recurrence cycle of a few years, a long duration 
of a few to several tens of days, a small amplitude up 
to 3 mag, a double-peaked structure in ultraviolet, the 
rise time equal to, or even longer than the decline time, 
and so on (Sabbadin, Bianchini 1983; Szkody et al. 1985; 
Bianchini ct al. 1986; Simon 2002; and references therein). 
Szkody et al. (1985) obtained optical spectra during an 
outburst (sec also Bianchini, Sabbadin 1981; Bianchini 
et al. 1982). He ii 4686 became stronger than in quies- 
cence, and the C iii/N ill Bowen blend appeared. The full 
widths at zero intensity (FWZI) of the emission lines were 
relatively small, ~40 A, almost equal to that in quiescence, 
which implied that the accretion disk did not extend down 
to the white dwarf surface either in quiescence or in out- 
burst. lUE ultraviolet spectra were acquired during the 
rise, at the maximum, and during the decline of the 1981 
outburst (Wu et al. 1989). These spectra were much red- 
der than observed in usual dwarf novae in outburst and in 
nova-likes, but Wu et al. (1989) noted that this is quali- 
tatively explained by a model of the accretion disk where 
the inner region is truncated by the magnetic field. Also 
in quiescence, black body fitting to an lUE spectra in- 
dicated that the temperature was not more than 12,000 
K, much lower than that of objects in this class (Rosino 
et al. 1982). Bianchini, Sabadin (1983) mentioned that 
the spectral energy distribution from X-ray to infrared of 
GK Per in quiescence is better explained by a model of a 
disk where magnetic fields control the temperature distri- 
bution, than another model of a standard disk with the 
inner region truncated. X-ray started brightening ^30 d 
earlier than optical light in the 1978 outburst (Bianchini, 
Sabbadin 1985). Based on this fact and that the UV out- 
burst light curve had two peaks, one of which occurred 
in the onset phase of the optical 1981 outburst, Bianchini 
et al. (1986) proposed that the outbursts in GK Per start 
at the inner edge of the accretion disk. 

A model to reproduce the peculiar outburst light curve 
was first proposed by Cannizzo, Kenyon (1986). They 
suggested, to explain the long outburst recurrence cycle 
and the long outburst duration, that a low viscosity pa- 
rameter a, about 1/10 of the typical value for dwarf novae 
with shorter Porb, and a very low mass transfer rate M 
of 10~^ AfQyr"^ were required. This condition makes an 
eruption due to the thermal instability occur at the inner 
edge of the accretion disk and the heating wave propagate 
to the outer disk. Such an outburst is called inside- out 
outburst. Although Angelini, Verbunt (1989) took into 
account the inner truncation of the accretion disk due 
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Fig. 3. Thcta diagrams and power spectra obtained using 
the combined data of February 27 and 28, March 3, 5, 6, 10, 
12, 13, 18, and 20, and April 5. (upper 2 panels) Periodgrams 
covering 0.25—1000 cycle d^^. There is no prominent peak in 
the range of >50 cycle d~^ . (lower 2 panel) Same as the upper 
panels, but those enlarged below 50 cycle d~^. Two peaks of 
0.05947(1) d (16.62 cycle d'^) and 0.05947(3) d (17.62 cycle 
d^^) have almost same significance. 

to the magnetic fields on the white dwarf, a low a was 
still required to produce the long duration of outburst. 
Kim et al. (1992) (see also Kim et al. 1992) proposed an- 
other model of the inside-out outburst requiring a larger 
inner radius (~2.5-4.5 x 10^" cm) of the disk and a large 
mass transfer rate (~ 10^^ g s"^), adopting a normal a 
value. Introducing the stagnation phenomenon enabled 
them to reproduce the double-peaked shape in UV ob- 
served by Bianchini et al. (1986). Yi et al. (1992), Yi, 
Vishniac (1994), and Yi, Kenyon (1997) developed this 
model to explain the spectral feature in the hard X-ray 
and ultraviolet range during optical outburst and dur- 
ing quiescence. Simon (2002) revealed strong relation- 
ship among the recurrence cycle, the outburst duration, 
and the outburst maximum magnitude, and that they had 
became longer, longer, and brighter, respectively, in the 
last five decades. These changes and that the final decay 
branch is very similar for all the outbursts may be a re- 
sult of three factors: decrease of irradiation of the disk by 
the primary star (cf. Schreiber, Gansicke 2001), decrease 
of viscosity, and consequently that the thermal instability 
starts at a different distance from the disk center for the 



D. Nogami et al. [Vol. 



-0.06 - 

-0.04 - 

0.02 - 

- 

a! - 

3 0.03 7 J 

S) 0.04 J * t 

S 1 1 — 

aj-0.06 - 

ld-0.04 - 

K 0.02 : 

'- - 

0.02 '- J 

0.04 i J ' 5 



P = 0.05613 d 



P = 0.05947 d 



i I I i ' I 



1 1- 



0.5 
Phase 



Fig. 4. Light curves folded with the periods of 0.05612 d and 
0.05947 d, using the data obtained on February 27 and 28, 
March 3, 5, 6, 10, 12, 13, 18, and 20, and April 5. 

respective outbursts. 

3. Observations 

We performed the observations at the Ouda Station, 
Kyoto University. A 60-cm reflector (focal length=4.8 m) 
and a CCD camera (Thomson TH 7882, 576 x 384 pix- 
els) attached to the Cassegrain focus were used (for more 
information of the instruments, see Ohtani et al. 1992). 
The on-chip 2 x 2, or occasionally 3x3, binning mode 
was selected to reduce the read-out and saving dead time. 
We adopted Johnson V- and i?-band interference filters. 
Table 1 gives the journal of the observation. The nominal 
error for each point is typically 0.01 mag. After standard 
de-biasing and flat fielding, the frames were processed by 
a microcomputer-based aperture photometry package de- 
veloped by one of the authors (TK). 

The magnitudes of the object were measured relative to 
a local standard star, HD 21588 {V^9M9, B -V^0.21Q, 
spectral type A2, in Hipparcos Input Catalogue. Version 
2 (Turon et al. 1993)). Heliocentric corrections to obser- 
vation times were applied before the following analysis. 

4. Results 

The long-term light curve drawn with visual observa- 
tions reported to VSNET^ and our observations around 
this outburst is exhibited in figure 1. The outburst is 
almost symmetric, but the rise time is somewhat longer 
than the decline time. The shape is divided into the fol- 
lowing three parts: 1) the slow rise with several times 
changes of the rising rate to the maximum for about 35 
d, 2) the subsequent gradual dechne with a rate of 20.0 d 
mag-i for about 16 d till HJD 2450188, and 3) the follow- 
ing rapid decline with a rate of 5.6 d mag~^ for about 10 
d. The outburst duration was about 60 d and the mag- 
nitude at the maximum is = 10.3. Although the 1996 

^ see <http;//www.kusastro.kyoto-u.ac.jp/vsnct/> 
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Fig. 5. Power spectra taken for the data on February 28, 
March 3, 5, and 6, to search for kilo-second QPOs. The QPO 
periods are measured to be ~5400 s (14.8 d"!), ~4300/'-3000 
s (20.2/28.8 d-i), ~4900 s (17.5 d"!), and 2800 s (30.6 d~^) 
on February 28, March 3, 5, and 6, respectively. 

outburst was not covered around the end because of the 
seasonal reason, a short fading tail probably followed the 
rapid decline phase (Simon 2002). 

GK Per became dramatically bluer for the first 5 days 
of the outburst, compared to its quiescence B — V color of 
0.9-1.0 (Bianchini et al. 1986; Wu et al. 1989; Patterson 
1991). The color reached its peak [B -V ^ 0.18) about 10 
d before the outburst maximum, then started to redden. 

Daily light curves are shown in figure 2. We can see 
modulations with amplitudes up to 0.5 mag. Selecting 
data of February 27 and 28, March 3, 5, 6, 10, 12, 13, 
18, and 20, and April 5 for their relatively long coverages 
(>1 hr); we performed a period analyses of the Phase 
Dispersion Minimization method (Stellingwcrf 1978) and 
the Fourier transfer, after subtraction of the daily linear 
trends, to detect oscillations existing for a long timescale. 
Figure 3 exhibits the resultant theta diagrams and power 
spectra. Two peaks indicating 0.05612(1) d (= 4849 s) 
and 0.05947(3) d (= 5138 s), which are one-day aliases of 
each other, have almost same significance. 

The fight curves folded with 0.05612 d and 0.05947 d, 
using the data same as those for the period analyses, are 
drawn in the lower and upper panels of figure 4, respec- 
tively. The full amplitudes are 0.10 mag, which is much 
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Fig. 6. Power spectra of daily pre-whitened data to seek short period modulations. The dotted lines indicate the spin frequency 
of 246 cycle d~^ (351 s). The lower arrows point at the frequency of 229 cycle d~^ (378 s), beat of which with the spin frequency 
gives rise to the mean QPO frequency of 5000 s. The spin feature is detected on March 13, and possibly on March 3. On March 3, 
and 6, QPOs with a typical period longer than the spin period were caught, while those with a typical period shorter than the spin 
period were detected on March 3, and 6, and April 5. 



smaller than that of apparent modulations in figure 2. The 
ephemerides are 

Ta{HJD) = 2450100.04336 + 0.05612 x E, (1) 

and 

Ta{HJD) = 2450100.03474 + 0.05947 x E. (2) 

The peaks calculated by equation (1) and (2) arc indicated 
by the tick marks and the upper arrows in figure 2. These 
marks do not necessarily correspond to observed maxima, 
which, in conjunction with the small amplitudes, means 
that the period of the modulations varies and the 0.05612- 
d or 0.05947-d period is just a mean period. Conclusively, 
QPOs with a typical period of ^5000 s were proved to 
exist during the current outburst. 

For closer watch of the evolution of the kilo-second 
QPOs, wc took power spectra of each onc-night data set 
covering over 3 hours. Figure 5 shows that those QPOs 
have a period of 5400 s, 4300/3000 s, 4900 s, and 2800 s, 
with a error of a few hundred seconds, on February 28, 
March 3, 5, 6, respectively. 



Wc again took power spectra of the each one- night data 
set of over-l-hour coverage, after pre-whitening to reveal 
short period modulations more clearly. The method of 
pre-whitening we adopted is to subtract, from each point, 
the magnitude calculated by averaging data within 7 min- 
utes before and after the point. The resultant power spec- 
tra are shown in figure 6. In the March 13th data, and less 
clearly in the March 3rd data we found, for the first time in 
optical photometry obtained during outburst, clear signs 
of oscillations with the spin period (350.5±2.0 s). We 
measured QPOs with a typical period of 413.0±1.6 s (209 
cycle d-i), and 456.9±1.6 s (189 cycle d'^) on March 
3, and 6, respectively, and marginally caught them with 
a period of 371.9±1.6 s (232 cycle d"^) on February 28. 
These periods arc longer than Pgpin- At the same time, 
QPOs with a typical period of 304.7±0.9 s (284cycle d'^), 
292.2±0.7 s (296 cycle d"!), and 290.3±2.6 s (297 cycle 
d~^) on March 3, and 6, and April 5, were revealed. Our 
observations first showed clear evidence of existence of os- 
cillations having a period shorter than the spin period. 
The profiles of these short QPOs were depicted in figure 
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Fig. 7. Short QPOs detected in figure 6 



7. 

The detected periods of short- and long-period QPOs 
are summarized in table 2 



5. Discussion 

5.1. short- and long-term QPOs 

As described in the previous section, we detected two 
periodicities for short-term QPOs, which are ~300s and 
'-370-460S, on February 28, March 3, 6, and April 5 dur- 
ing our observation since 1996 February 26. Since the 
amplitude of the QPOs varies with a timescale of hours 
(Patterson 1991; Mazeh et al. 1985), it is hkely that QPOs 
had amplitudes smaller than the detection limit on the 
other nights. 

Morales-Rueda et al. (1999) spectroscopically observed 
GK Per on 1996 February 26, 27, and 28, and reported re- 
sults of detailed period analyses. While kilo-second QPOs 
and spin pulsations were detected in the V/R ratio vari- 
ation of H/3 and He ii 4686 in all nights except for spin 
pulsations in He ii on February 26, they also caught kilo- 
second QPOs in the integrated line fluxes of H/? and He 
II and the continuum in all nights. Spin pulsations in 
the integrated line fluxes were found only on February 



26. The typical period of the kilo-second QPOs of the 
continuum around H/? and He ii 4686 in Morales-Rueda 
et al. (1999) was 4477 (± 12) s (= 0.0518 d; 19.3 cycle d'^) 
on February 27 and 3692 (± 9) s (= 0.0427 d; 23.4 cycle 
d^^) on February 28. However, spin pulsations were not 
detected in our data on February 27, and 28. Contribution 
of the continuum would dilute the spin pulsations. The 
5400-s period we detected on February 28 was far longer 
than 3692 s found by Morales-Rueda et al. (1999) on the 
same day. This may suggest that the QPO period depends 
on the wavelength. 

The period of the kilo-second QPOs observed by 
Morales-Rueda et al. (1999) decreased from night to night 
in February 26-28, and they interpreted that this might 
be due to shrinkage of the blobs' orbit. The kilo-second 
QPO period, however, once increased on March 5 (table 
2). This period increase may be explained by generation 
of a new blobs on the outer orbit, but note that the period 
rapidly decreased to ~2800 s by the next night. 

Concerning short-term QPOs, we found two typical pe- 
riods of ^300 s and ~440 s in V band. The current 
standard scheme of QPOs in GK Per is the beat fre- 
quency model which was originally proposed by Lamb 
et al. (1985) for QPOs in the X-ray binary Sco X-1. 
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For QPOs in GK Per, Watson et al. (1985) first inter- 
preted the X-ray QPOs with this model that short-term 
QPOs are originated in dense blobs orbiting near the inner 
edge of the accretion disk with the Keplerian frequency 
and long-term QPOs observed in X-ray and optical has 
the beat frequency between the spin frequency and the 
blobs' frequency. Developing this model, Morales-Rueda 
et al. (1999) proposed a disk-curtain beat mechanism. The 
other plausible model is the disk-overflow accretion model 
proposed by Hellier, Livio (1994). This model is based on 
the idea that the kilo-second period is the characteristic 
timcscale of the radius where an accretion stream over- 
flowing the disk would collide back onto the disk. 

Simple calculation by the beat model yields the long- 
term QPO period of 1,500 s-2,400 s for both of -440- 
s and — 300-s periodicities we observed. These periods, 
however, can not be reconciled with the period of 2800- 
5400 s in our data, which might favor the disk-overflow 
accretion model. But the coverages of our data were not 
neccesarily long enough for accurate determination of the 
kilo-second periods. It is important to organize a multi- 
longitudinal observation campaign to simultaneously in- 
vestigate periodicities of hundreds of seconds and kilo- 
seconds. Note that using a bluer band is better in this 
case, since the amplitude of QPOs is larger in a bluer 
band (see Patterson 1991). 

Oscillations with a period of —300 s are the first dis- 
covery of below-the-Pspin modulations except for the 40-s 
pulse reported by Pezzuto et al. (1996) whose origin is 
still unknown. If this period is the Keplerian period of 
blobs, the inner edge of the accretion disk approached the 
white dwarf surface over the co-rotation radius during the 
1996 outburst, while the inner radius in quiescence is re- 
quired to be far larger than the co-rotation radius (Kim 
et al. 1992). The QPO period in fact decreased from 305 
s to 290 s in March 3-April 5. Theoretical exploration 
on how closely the inner radius can approach the primary 
star against the magnetic torque during outburst should 
be an interesting topic. It should be also pointed out that 
two periodicities were caught on the same nights: 305 s 
and 413 s on March 3 and 457 s and 292 s on March 6. 

5.2. Outburst mechanism 

Outbursts in dwarf novae are currently thought to be 
caused by the disk instability. However, the outburst 



mechanism in GK Per is still an open problem, because 
of the IP nature, the peculiar binary parameters such 
as the quite long Porb, and the extraordinary outburst 
property shown in figure 1 and summarized in table 3. 
While Bianchini et al. (1982) explained the 1981 out- 
burst by an increase in the mass transfer rate from the 
late-type secondary, Cannizzo, Kenyon (1986), Angelini, 
Verbunt (1989), and Kim et al. (1992) proposed models 
based on the disk instability scheme, as mentioned in sec- 
tion 2. 

SS Cyg-type dwarf novae usually show asymmetric out- 
bursts with a rapid rise and a slow decay, but some stars 
occasionally cause rather symmetric outbursts with slow 
rise, say, "anomalous" outbursts (for example, Cannizzo, 
Mattel 1992 and references therein for the SS Cyg case). 
Theoretical investigations indicate that the former is the 
outburst of the outside-in type and the anomalous out- 
burst is the one of the inside-out type (Smak 1984). The 
outbursts in GK Per are very similar to the anomalous 
outburst in shape, although the timescales of rise and 
decline are much longer in GK Per than in usual dwarf 
novae. Therefore, Cannizzo, Kenyon (1986), Angelini, 
Verbunt (1989), and Kim et al. (1992) attempted to apply 
an idea of the inside- out outburst for explanation of the 
outburst pattern of GK Per. 

Our observations revealed that the B — V color reached 
its minimum about 10 days before the outburst maxi- 
mum. This agrees with the fact that SS Cyg was bluest 
in the B — V color 4.0 d before the F-band maximum 
during the 1981 September anomalous outburst (Hopp, 
Wolk 1984), though SS Cyg become bluest at, or a lit- 
tle later than the maximum during normal outbursts (e.g. 
Hopp, Witzigmann 1980). The B — V observations pre- 
sented here support that outbursts in GK Per arc of the 
inside-out type. 

Although Cannizzo, Kenyon (1986) did not take into ac- 
count the inner-disk truncation due to the magnetic fields 
on the white dwarf, their calculation provided a prediction 
that GK Per becomes bluer in B ~V for a while after the 
outburst maximum, which clearly contradicts with our re- 
sult. In polishing up the model for GK Per being rooted 
in the disk instability theory, it may be a key whether the 
model can account for the B — V variation during out- 
burst. 

Kim ct al. (1992) proposed a time-dependent disk insta- 



10 



D. Nogami ct al. 

Table 3. Properties of IPs showing outbursts. 
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Note: In some case, we derived the values of outburst duration, outburst amplitude, and the recurrence cycle from 
the observations reported to VSNET. 

* In only one case, an outburst was caught just 8 days after the previous one. 

^ These large amplitude outbursts occur only when TV Col is in high state (Augusteijn et al. 1994). 
Reference: 1 Tovmassian et al. (1998), 2 Ishioka et al. (2002), 3 Helher, Sproats (1992), 4 HeUier et al. (2000), 
5 van Amcrongen et al. (1987), 6 Jablonski, Steiner (1997), 7 Wenzel (1983), 8 Mateo et al. (1991), 9 Patterson, 
Szkody (1993), 10 HaswcU ct al. (1997), 11 Norton et al. (1999), 12 Simon (2000), 13 Szkody et al. (2002), 14 Schrijvcr 
et al. (1985), 15 HcUicr, Buckley (1993), 16 Augusteijn ct al. (1994), 17 Takano ct al. (1989), 18 Kamata et al. (1991), 
19 Hellier ct al. (1997), 20 Watson et al. (1985), 21 Morales-Rucda ct al. (2002), 22 Simon (2002) 



bility model to reproduce the recurrence time, the dura- 
tion, the doubly-peaked structure in the UV range. This 
model appears to be able to explain also the outburst 
shape consisting of the three parts (see e.g. figure 6 in 
their paper), although the radius of the inner edge must 
be changed for the individual outburst in order to repro- 
duce the variation of the rise-branch shape. Simon (2002) 
found that the rise branch had showed large scatter while 
the decline rate during the final rapid decline was common 
to 8 outbursts in these three decades. 

We here focus on the final rapid decline. During this 
part, the cooling front from the outer region propagates 
onto the inner edge (see figures in Kim et al. 1992). This 
process is common to the outbursts of the outsidc-in type 
and of the inside-out type, so that the decline rate dur- 
ing this phase is expected to be on the Bailey's relation 
(Bailey 1975; Szkody, Mattel 1984; Warner 1995). The 
value of 5.6 d mag~^ in figure 1, however, is much smaller 
than 10.0 d mag~^ calculated using the Bailey's relation 
for the Porb of GK Per. Simon (2002) derived 8.0 d mag^^ 
from the recent 8 outbursts [they also derived a value 
of 14.5 d mag~^ using the data including those during 
the gradual fading. However, this treatment is not ad- 
equate because different processes are thought to work 
during the gradual fading phase and during the rapid fad- 
ing phase]. Also in other IPs, DO Dra (Simon 2002) and 
HT Cam (Ishioka et al. 2002), the decay timescale is re- 
ported to be shorter than expected from their Porb by the 
Bailey's relation. The decay timescale is made shorter by 
the inner truncation of the disk (Angclini, Verbunt 1989; 
Cannizzo 1994), and the magnetic fields penetrating the 
accretion disk may also have an effect on the timescale. 

Properties of GK Per and other outbursting IPs are 
summarized in table 3. We can not see a simple corre- 
lation between any two factors among 5 systems. This 
means that it is difficult to establish a unified scheme on 



the outburst mechanism in IPs. Note that short outbursts 
in V1223 Sgr and TV Col are possibly caused by the mass 
transfer events, not by the disk instability (e.g. Hellier, 
Buckley 1993). A modeling including the effects of mag- 
netic fields penetrating and truncating the accretion disk, 
3D structures of the accretion disk, the mass stream, the 
accretion curtain, and the accretion pole, will be required 
to reproduce a variety of outbursts in IPs. 

6. Conclusion 

Our I^-band timc-rcsolved photometry revealed the spin 
pulse (Pspin = 351 s), and QPOs with three typical 
timescales of 300 s, 440 s, and 5,000 s. The discovery 
of 300-s QPOs means that the inner edge of the accretion 
disk shrinks over the co-rotation period during outburst. 
The beat period of 440 s and 351 s is —1,700 s, which 
would be against the model for the kilo-second QPOs. 
Further check of contemporaneousness of hundred-second 
and kilo-second QPOs is required. The B — V color be- 
came bluest (P — 1/ — 0.18) about 10 d before the outburst 
maximum. A similar color variation was observed during 
an anomalous outburst in SS Cyg, and this agrees with the 
idea that GK Per shows inside-out- type outbursts. The 
decline rate during the rapid decline phase is smaller in 
GK Per, as well as DO Dra, and HT Cam, than expected 
from the Porb- This is qualitatively explained by trun- 
cation of the inner part of the accretion disk, but more 
detailed modeling taking into account the effect of the 
magnetic fields, 3D structures of the accretion disk, the 
mass stream, the accretion curtain, and the accretion pole 
is needed to reproduce the whole outburst properties of 
IPs. 

We thank the anonymous referee for the useful com- 
ments. The authors are very grateful to amateur observers 
for continuous reporting their valuable observations to 
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of Education, Culture, Sports, Science and Technology 
(TK). 
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